Measurement of da/dM and Forward-Backward Charge Asymmetry for High Mass 
Drell-Yan e + e Pairs from pp Collisions at yfs = 1.8 TeV 
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We report on a measurement of the mass dependence of the forward-backward charge asymmetry, 
Afb, and production cross section da/dM for e e~ pairs with mass M ee > 40 GeV/c 2 . The data 
sample consists of 108 pb _1 of pp collisions at ^/s = 1.8 TeV taken by the Collider Detector at Fer- 
milab during 1992-1995. The measured asymmetry and da/dM are compared with the predictions 
of the Standard Model and a model with an extra Z' gauge boson. 

PACS numbers: 13.85.Qk, 12.38.Qk, 12.15.Ji, 12.15.Mm, 12.60.-i, 12.60.Cn 



In hadron-hadron collisions at high energies, massive 
e + e~ pairs are predominantly produced via the Drell- 
Yan process. In the Standard Model, quark-antiquark 
annihilations form an intermediate 7* or Z (7* /Z) vec- 
tor boson, which may yield an e + e" pair. The presence 
of both vector and axial- vector couplings in this process 
gives rise to a forward-backward asymmetry, Afb, in the 
final-state angle of the electron in the rest frame of the 
e + e _ pair (with respect to the proton direction). The 
Standard Model predicts accurately the mass [M] depen- 
dence of the production cross section (da/dM) and Afb- 
In the region of the Z resonance, the predicted asym- 
metry is related to the electroweak mixing angle. For 
M ^> Mz, the large predicted asymmetry (near 0.6) is 
a consequence of the interference between the propaga- 
tors of the 7* and Z. Various new interactions not in- 
cluded in the Standard Model could result in deviations 
from the Standard Model predictions for both da/dM 
and Afb- Possible new interactions include additional 
gauge bosons |@-||], quark- lepton compositeness M, ex- 
change of R parity violating supersymmetric particles or 
leptoquarks, and extra-dimensions ||. 

In this communication we extract measurements of 
da/dM and improve the existing measurements of Afb- 
The previous CDF measurements in the e + e~ [gj and 
@ channels reported on a central (\rj\ < 1.1) re- 
gion measurements of d 2 a/dMdy averaged over rapidi- 
ties H of \y\ < 1. Electrons in the plug (1.1 < \rj\ < 2.4), 
but not forward (2.2 < |^| < 4.2) calorimeter were in- 
cluded in the previous CDF measurement || of Afb, 
but backgrounds in the plug region were large resulting 
in significant uncertainties. This analysis includes elec- 
trons in the plug and forward regions of the calorimeter, 
thus covering a larger range of rapidity. New analysis 
techniques and tracking requirements in the forward re- 
gion greatly reduce the backgrounds and uncertainties. 
In addition, by applying electroweak and QED radiative 
corrections [ |10|Jir| (which correct for the change in M 
from the emission of final state photons), measurements 
of Afb in small bins over a large range in mass are ex- 
tracted for the first time. Our measurement of e + e~ 
pairs with invariant masses between 40 and 500 GeV/c 2 
complement measurements of cross-sections and asym- 
metries in e + e~ machines, which currently extend to 200 
GeV/c 2 . 

The e + e~ pairs come from 108 pb _1 of pp collisions 
at y/s = 1.8 TeV taken by the Collider Detector at Fer- 
milab § (CDF) during 1992-1993 (18.7 ±0.7 pb" 1 ) and 



1994-1995 (89.1 ± 3.7 pb" 1 ). The CDF detector con- 
sists of a solenoidal magnetic spectrometer surrounded by 
projective-tower-geometry calorimeters and outer muon 
detectors. Only detector components used in this mea- 
surement are described here. Charged particle momenta 
and directions are measured by the spectrometer embed- 
ded in a 1.4 T axial magnetic field, an 84-layer cylindrical 
drift chamber (CTC), an inner vertex tracking chamber 
(VTX), and a silicon vertex detector (SVX). The polar 
coverage of the CTC tracking is \rj\ < 1.2. The pp col- 
lision point along the beam line (Z veTtex ) is determined 
using tracks in the VTX. The energies and directions || 
of electrons, photons, and jets are measured by three sep- 
arate calorimeters covering three regions: central, end 
plug, and forward. Each calorimeter has an electromag- 
netic (EM) and a hadronic (HAD) calorimeter. We fol- 
low the analysis procedure used in our previous measure- 
ments fl2] , |i~3|| of da/dy and da/dP^ of Z boson pairs. 
The measurement of Afb requires the sign of at least 
one lepton. Therefore, we require one of the leptons to 
be in the central region. 

The sample of e + e~ events was collected by a three- 
level online trigger that required an electron in either the 
central or the plug calorimeter. The offline analysis se- 
lected events with two or more electron candidates. Since 
the electrons from the Drell-Yan process are typically iso- 
lated, a cut on the amount of transverse energy in a cone 
outside the electron shower is imposed. Electrons in the 
central, end plug, and forward regions are required to be 
within the fiducial area of the calorimeters and to have 
a minimum Et of 22, 20, and 15 GeV, respectively. To 
improve the purity of the sample, electron identification 
cuts are applied H] . The central electron (or one of them 
if there are two) is required to pass strict criteria. The cri- 
teria on the other electron are looser. A central electron 
must have a CTC track that extrapolates to the elec- 
tron's shower cluster in the EM calorimeter. This cluster 
must have EM-like transverse shower profiles. The track 
momentum and the EM shower energy must be consis- 
tent with one another. The track is also used to deter- 
mine the position and direction of the central electron. 
The fraction of energy in the HAD calorimeter towers 
behind the EM shower is required to be consistent with 
that expected for an EM shower (Shad/^em)- The plug 
electrons must also have an EM-like transverse shower 
profile. The end plug and forward electrons are required 
to pass the -Ehad / Eem requirement and to have a track 
in the VTX which originates from the same vertex as the 
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other electron and points to the position of the electro- 
magnetic cluster in the calorimeter. The ratio of found 
to expected hits in the VTX is required to be greater 
than 70% and 50% for plug and forward electrons, re- 
spectively. The VTX tracking efficiency is (97.8 ± 0.3)% 
for plug electrons and (97.0±0.9)% for forward electrons. 
To further reduce the backgrounds from dijet events for 
very forward electrons both leptons legs for CF events 
are required to be on the same side of the detector i.e. 
have the same sign in r\. The total data sample after all 
cuts consists of 7632 events. 

With the above cuts, backgrounds are low and can be 
reliably estimated using the data. Because of the CTC 
tracking requirement, the jet background for the CC sam- 
ple is negligible. All CC events are required to be oppo- 
site sign (the 8 same sign events are all in the Z mass 
region and originate from the small charge misidentifi- 
cation probability). The CC background is mainly com- 



cc, bb, and tt 



posed of e + e~ pairs from W + W~ 

sources. This background is estimated using e ± ^ pairs 
Because of the tracking and same side requirements, the 
backgrounds in the CP and CF topologies are also small 
and mainly come from dijets events. These are extracted 
from the data by fitting the isolation energy distribu- 
tion to a sum of the expected distributions for signal 
and background events. For M ee > 105 GeV the isola- 
tion energy fraction is used. The shapes of the isolation 
energy distributions are extracted from specially selected 
jet background-dominated and electron signal-dominated 
data samples. 

The acceptance for Drell-Yan e + e~ pairs is ob- 
tained using the Monte Carlo event generator PYTHIA 
6.146 and CDF detector simulation programs. 

PYTHIA generates the LO QCD interaction (q + q — > 
7* /Z), simulates initial state QCD radiation via its par- 
ton shower algorithms, and generates the decay, 7* jZ — > 
e + e~. To approximate higher order QCD corrections 
to the LO mass distribution, a X-factor pi is used as 
an event weight: K(M 2 ) = 1 + f (1 + ^TT*)a s (M 2 )/2TT : 
where a s is the QCD coupling. This factor improves the 
agreement between the NLO and LO mass spectra. (For 
M > 50 GeV/c 2 , 1.25 < K < 1.36.) The CTEQ5L @ 
PDFs are used in the acceptance calculations. Final state 
QED radiation |l(J from the 7* /Z — > e + e~ vertex is 
added by the PHOTOS |l|] Monte Carlo program. The 
parameters of PYTHIA 6.146 are tuned Il5^to fit experi- 
mentally measured transverse momentum] 13j and rapid- 
ity distributions |l2]| of e + e~ pairs in the Z region. Gen- 
erated events are processed by CDF detector simulation 
programs and are reconstructed as data. The calorimetry 
energy scales, resolutions, cut efficiencies and correspond- 
ing errors used in the detector simulation are extracted 
from the data. Reconstructed Monte Carlo events are 
accepted if they pass the fiducial, kinematic, and mass 
cuts. 

The angle of the electron is determined in the Collins- 
Soper jl8| frame, which reduces the uncertainty in the an- 
gle introduced by the transverse momentum of the incom- 



ing partons. We determine the cross sections da + /dM 
and da^ /dM for forward and backward events, respec- 
tively, da / dM (shown in Figure 0(a)) and Afb (shown in 
Figure 0(b)) are extracted from the sum and difference of 
the forward and backward cross sections. The very small 
misidentification probability is handled as a dilution fac- 
tor in Afb- The CC, CP, and CF samples are combined 
and are binned in M, and da j dM is given by 



da 
dM 



N - B 



C rc AAd y ] r £ r cA r 



Here, N is the number of observed events, B is the esti- 
mated background in a bin, C rc is the radiative correc- 
tion |H],[ll| factor to convert the measured cross sections 
to a Born level (one boson exchange) cross section, AAf 
is the bin width, the sum r is over the 1992-1993 and 
1994-1995 runs, C r is the integrated luminosity, and eA r 
is the run's combined event selection efficiency and ac- 
ceptance. The backgrounds subtracted from the event 
count are predicted using the data and background sam- 
ples. Acceptances are calculated separately for CC, CP, 
and CF pairs. They are combined with the correspond- 
ing event selection efficiencies to give eA r . A r includes 
corrections for detector resolution smearing. Since the 
PYTHIA/PHOTOS (H) Monte Carlo includes radiative 
final state photons, the effects of the final state radi- 
ated photons are calculated together with the accep- 
tance calculation. Also shown in Figure |l](a) are CDF 
dimuon results for da / dM (which have not been previ- 
ously published). These were extracted from the pub- 
lished dimuon Q data (d 2 a/dMdy for \y\ < 1) byapply- 
ing a correction using a QCD-NLO calculation [[19| that 
describes our measured Z e + e~ pair y distributions |l2j ] 
over the entire kinematic range. The e + e~ cross sections 
and asymmetries for M ee < 105 GeV/c 2 are given in Ta- 
ble ||, and the higher mass data samples and backgrounds 
are summarized in Table 0. We define M to be the mass 
value for which the cross section is equal to the average 
cross section over each bin, as calculated from a NNLO 
theoretical prediction p(J] for da/dM. 

The systematic errors are determined from the follow- 
ing: statistical errors in the fits to the fractional back- 
ground, variations in the background estimates using dif- 
ferent methods, the background in the efficiency sample, 
the uncertainty in energy resolution of the calorimeter, 
the choice of PDFs, and the distribution of e + e~ Pt used 
in the Monte Carlo event generator. Above the region of 
the Z mass, the total systematic errors are much smaller 
than the statistical errors, and the total errors shown 
in Table || are dominated by statistics. In the Z mass 
region, the systematic errors are less than half of the 
statistical error (as shown below). The pp collision lumi- 
nosity is derived with CDF's beam-beam cross section, 
cr BBC = 51.15 ± 1.60 mb The luminosity error of 

3.9% contains the itbbc error and uncertainties specific 
to running conditions. As described in a previous com- 
munication [Ol , the extracted cross section in the Z mass 
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region is 252. l±3.9(stat.)±1.6(syst.)±9.8(lum.) pb. Note 
that since the pp inelastic cross section used by CDF in 
luminosity calculations differs from D0's by +5.9% [f22| 



the DO e + e cross sections |2^] shown in Figure [j] have 
been multiplied by 1.059. 
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FIG. 1. (a) da/dM distribution of e + e~ and pairs. 
All errors (except for the overall 3.9 % luminosity error) have 
been combined in quadrature. The Standard Model theoreti- 
cal predictions (solid line) have been normalized by a factor 
of 1.11 to the data in the Z boson mass region. Also shown 
are e + e~ measurements from D0. (b) Afb versus mass com- 
pared to the Standard Model expectation (solid line). Also, 
predicted theoretical curves for da/dM and Afb with an ex- 
tra E 6 Z' boson (width of 10%) with M z ,=350 (dotted line) 
and 500 GeV/c 2 (dashed line). The inset in (a) shows the 
difference between the da/dM CDF data for e + e~ pairs and 
the Standard Model prediction (on a linear scale) compared 
to expectation from these two Z' models. 

Figure [j] also compares the measured da/dM and Afb 
to theoretical predictions. The da/dM curve is a QCD 
NNLO fH calculation with MRST99 NLO PDFs ||. 
The predictions in Figure |l|(a) are normalized by a factor 
F=l.ll, the ratio of measured total cross section in the Z 
region to the NNLO prediction (the overall normalization 
uncertainties are 3.9% for the experimental data and 5% 
for the NNLO theory). The Standard Model predictions 
for A FB have been calculated §[l0) in NLO-QCD. The 
measured da/dM and Afb values are in good agreement 
with the Standard Model predictions. However, Afb in 
the highest mass bin (300-600 GeV/c 2 ) is 2.2cr (stan- 



dard deviations) below the Standard Model prediction. 
There are three events pj| in the negative hemisphere 
and one event in the positive hemisphere. A negative 
asymmetry in this region could result from new interac- 
tions including additional gauge bosons |^|-[| (discussed 
as an example below), quark-lepton compositeness [^), 
exchange of R parity violating supersymmetric particles, 
leptoquarks, and extra-dimensions ||. For example, the 
predicted theoretical curves j|] for da/dM and Afb in 
models which include extra E6 Z' bosons (with parame- 
ters tuned to fit low energy electroweak data) are shown 
as the dashed and dotted curves in Figure |l|. 

In summary, the measured da / dM and Afb values are 
in good agreement with the Standard Model predictions. 
At the highest mass values, the measured value of Afb 
is 2 a lower than the Standard Model prediction. These 
high mass Drell-Yan data will be included in global fits to 
electroweak data to search (or extract limits) for physics 
beyond the Standard Model. 
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TABLE I. The e + e~ cross sections and asymmetries for 
M ee < 105 GeV/c 2 . The statistical and systematic errors 
have been combined in quadrature. The 3.9% overall lumi- 
nosity error is not included. 



Mass Bin 


M ee 


da ee /dM 
(pb/GeV/c 2 ) 


A FB 


(GeV/c 2 ) 


(GeV/c 2 ) 




40-50 


44.5 


2.30 ±0.47 


-0.02 ±0.17 


50-60 


54.6 


0.80 ±0.11 


-0.17 ±0.10 


60-70 


63.8 


0.60 ±0.06 


-0.56 ±0.07 


70-78 


74.7 


0.58 ±0.05 


-0.47 ±0.06 


78-86 


82.9 


1.64 ±0.06 


-0.21 ± 0.04 


86-88 


87.1 


5.89 ±0.22 


-0.09 ±0.04 


88-90 


89.2 


18.28 ±0.52 


-0.04 ±0.03 


90-92 


90.6 


54.85 ± 1.39 


0.07 ±0.02 


92 - 94 


92.8 


22.75 ±0.66 


0.11 ±0.03 


94 - 100 


96.2 


3.80 ±0.12 


0.28 ±0.03 


100 - 105 


102.2 


0.91 ±0.07 


0.38 ± 0.07 



TABLE II. The number of signal and background events 
and correction factors involving the extraction of da^ e /dM, 
da^e/dM, daee/dM, and da^^/dM for mass bins above the Z 
pole. Here, N is the number of events, B is the background 
estimate, eAC rc is the combined efficiency, acceptance, smear- 
ing and radiative correction. The statistical and systematic 
errors for da/dM have been combined in quadrature. The 
3.9% overall luminosity error is not included. 



Mass Bin 


Mee 


N 


B 


eACrc 


date/dM 
(fb/GeV/c 2 ) 


(GeV/c 2 ) 


(GeV/c 2 ) 








105- 


120 


111.2 


93 


2.3 


0.293 


191.55 ±20.38 


120- 


140 


128.8 


34 


1.4 


0.307 


49.20 ± 8.81 


140- 


200 


164.2 


35 


1.2 


0.326 


-i r* f\~\ i c\ c\r\ 

16.01 ± 2.80 


200- 


300 


240.6 


7 


0.1 


0.337 


1.88 ± 0.73 


300 - 


400 


343.2 


1 


0.0 


0.330 


0.28 ± 0.28 


400 - 


600 


478.8 





0.0 


0.307 


0.00 ± 0.15 


600- 


999 


725.6 





0.0 


0.227 


0.00 ± 0.10 


Mass Bin 


Mee 


N 


B 




da^e/dM 


105- 


120 


111.2 


39 


2.3 


0.321 


70.64 ± 12.04 


120- 


140 


128.8 


13 


1.4 


0.334 


16.07 ± 5.01 


140- 


200 


164.2 


12 


1.2 


0.353 


4.74 ± 1.52 


200- 


300 


240.6 


3 


0.1 


0.386 


0.69 ± 0.42 


300 - 


400 


343.2 


2 


0.0 


0.398 


0.47 ± 0.33 


400 - 


600 


478.8 


1 


0.0 


0.391 


0.12 ± 0.12 


600 - 


999 


725.6 





0.0 


0.355 


0.00 ± 0.07 


Mass Bin 


Mee 


N 


B 




daee/dM 


105- 


120 


111.2 


132 


4.7 


0.300 


262.29 ± 23.67 


120- 


140 


128.8 


47 


2.9 


0.313 


65.46 ± 10.17 


140- 


200 


164.2 


47 


2.4 


0.331 


20.81 ±3.20 


200- 


300 


240.6 


10 


0.3 


0.346 


2.60 ±0.85 


300 - 


400 


342.2 


3 


0.0 


0.343 


0.81 ±0.47 


400 - 


600 


478.8 


1 


0.0 


0.324 


0.14 ±0.14 


600- 


999 


725.6 





0.0 


0.252 


0.00 ±0.09 


Mass Bin 


M w 


N 


B 




da^/dM 


110- 


120 


114.5 


29 





0.181 


231.93 ±51.56 


120- 


150 


132.5 


28 





0.167 


80.61 ± 15.93 


150- 


200 


171.2 


9 





0.164 


13.98 ± 4.67 


200- 


300 


240.6 


6 





0.168 


3.91 ± 1.58 


300- 


400 


343.2 


1 





0.176 


0.56 ±0.56 


400- 


600 


478.8 








0.227 


0.00 ±0.20 


600- 


999 


725.6 








0.230 


0.00 ±0.11 



7 



